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Abstract 

We consider in mean-field theory the continuous development below a second- 
order phase transition of n-atic tangent plane order on a deformable surface 
of genus zero with order parameter ip = [e tnd ). Tangent plane order expels 
Gaussian curvature. In addition, the total vorticity of orientational order on 
a surface of genus zero is two. Thus, the ordered phase of an n-atic on such 
a surface will have 2n vortices of strength 1/n, 2n zeros in its order param- 
eter, and a nonspherical equilibrium shape. Our calculations are based on 
a phenomenological model with a gauge-like coupling between tp and curva- 
ture, and our analysis follows closely the Abrikosov treatment of a type II 

superconductor just below H C 2- 
PACS numbers: 87.20E, 02.40, 64.70M 
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When aliphatic molecules are dissolved in water, they can spontaneously segregate into 
bilayer membranes in which hydrocarbon tails are shielded from contact with water. De- 
pending on conditions, these membranes can form random, extended surfaces, regular pe- 
riodic structures, or closed vesicles separating an interior from an exterior [1]. Membranes 
of a similar nature but with more structure and greater complexity form biological cell 
walls. In a wide class of membranes molecules move freely within the membrane forming 
a 2-dimensional fluid offering little resistance to changes in membrane shape. Such mem- 
branes are physical examples of random surfaces, which can undergo violent shape changes. 
Molecules in membranes can also exhibit varying degrees of orientational and positional 
order[2] including tilt, hexatic, and crystalline order similar to that found in free standing 
liquid crystal films [3]. These ordered membranes provide fascinating laboratories for the 
study of the coupling between order and geometry, analogous in many ways to the coupling 
between matter and geometry in general relativity. The underlying cause of this coupling is 
easy to see. A vector field (or any field describing orientational order) cannot be everywhere 
parallel if it is forced to lie on a surface, such as that of a sphere, that curves in two directions 
and thus has nonzero Gaussian curvature. It could lower its energy by flattening the surface. 
An additional complication arises when order develops on closed surfaces. A closed surface 
can be classified according to its genus (or number of handles) g: a sphere has genus zero, a 
torus genus one, etc. Orientational order on a closed surface necessarily [4] has topological 
defects (vortices) with total strength (vorticity) equal to the Euler characteristic 2(1 — g) of 
the surface. Tangent plane order on a sphere will have vorticity 2, a torus vorticity 0, etc. 
The continuous development of vector order on a deformable surface of genus zero will be 
accompanied by a continuous change from spherical to ellipsoidal shape[5]. Since vortices 
are energetically costly, it may be favorable for a closed physical membrane to transform into 
an open cylindrical structure [6] when tangent plane order develops in response to changes in 
temperature or other control variable. Indeed, there are a number of experimental examples 
of shapes changes [7] that may be explained by the development of tangent plane order. 

In this paper, we investigate in mean-field theory the development of n-atic order on a 
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closed surface of genus zero and the concomitant change in shape from spherical to non- 
spherical. An n-atic order parameter can have vortices of strength 1/n, and, since it is 
generally favorable to form vortices of minimum strength, we expect 2n maximally separated 
vortices of strength 1/n to be present in the ordered phase. Thus for n — 1, we expect two 
antipodal vortices[5], and respectively for n = 2, n = 3 and n = 6, we expect vortices at 
the vertices of a tetrahedron, an octahedron, and an icosahedron. Indeed calculations on a 
rigid sphere confirm this conjecture [6, 8]. We also expect the vesicle shape to change from 
spherical to ellipsoidal, tetrahedral, octahedral, or icosahedral in the four cases above. Fig. 
1 shows our calculated shapes for temperatures just below the transition temperature to the 
n-atic state for n = 1, 2, 3, 4, and 6. 

Our calculations are based on a phenomeno logical Hamiltonian[9] for a complex order 
parameter field whose coupling to shape occurs via a covariant derivative and via changes 
in the metric tensor. The model is almost identical to the Landau-Ginzburg theory of 
superconductivity except that vorticity is fixed by surface topology rather than energetically 
determined by an external magnetic field. The ordering transition we find is very similar to 
the transition from a normal metal to the Abrikosov vortex lattice [10] in a superconductor, 
and indeed our analysis follows very closely that of Abrikosov. We find the order parameter, 
which is necessarily spatially inhomogeneous because of the topological constraint on the 
total vorticity, from among a highly degenerate set of functions that diagonalize the harmonic 
Hamiltonian on a rigid sphere. This degenerate set has exactly 2n zeros at arbitrary positions 
on the sphere and is very similar in form to fractional quantum Hall wave functions [11]. 

Positions on a two-dimensional surface in R 3 are specified by a vector R(cr) as a function 
of a two-dimensional coordinate cr = (o"i,o- 2 ). Associated with R(<x) is a metric tensor 
Qabi. ') = d a H((r)-dbR((T) and a curvature tensor K ab (cr) defined via d a dbR = K ab N where N 
is the local unit normal to the surface. To describe tangent-plane n-atic order, we introduce 
orthonormal unit vectors ei and e 2 at each point on the surface. If t(<x) is any unit tangent 
vector, then e^cr) -t(cr) = cos#(cr) defines a local angle 9(a). n-atic order is then described 
by the order parameter ip(cr) = ^ e m6l ( cr )^ which can be related to the nth rank symmetric 
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traceless tensor constructed from the vector t. Note that since 9(a) depends on the choice of 
orthonormal vectors e x and e 2 , the order parameter ip(cr) does as well. This means that any 
spatial derivatives in a phenomenological Hamiltonian for ip must be covariant derivatives. 
The long-wavelength Landau- Ginzburg Hamiltonian for ip is Ti^ = Hi + H2 where 



where A a = e 1 ■ d a e 2 is the spin connection, g = detg ab and g ab gb c = When n = 6, this 
H^, is identical to the Hamiltonian for hexatic order on a deformable surface introduced in 
Ref. [12]. It is important to note that the covariant derivative D a = d a — inA a depends on 
n. The energy associated with shape changes of a constant area membrane is described by 
the Helfrich Hamiltonian [13], 



The complete Hamiltonian for n-atic order on a deformable surface is H = + 7i CU rv 

We can now specialize to surfaces of genus zero with constant area A = AirR 2 whose 
shapes do not differ significantly from that of a sphere. We choose cr = (6, 0) = Q to be 
the polar coordinates of a sphere and set R(f2) — R [l + p(f2)]e r , where e r is the radial unit 
vector. This allows us to introduce reduced metric and curvature tensors g ab and K ab via 
g a b = Rodab an d Kab = R^^ab- The field p(fi) measures deviations from sphericity and can 
be expanded in spherical harmonics. Any isotropic change in R can be described by R . In 
addition, uniform translations, which change neither the shape or the energy of the vesicle, 
correspond to distortions in p with / = 1. These considerations imply that p will have no 
I — or I — 1 components: p(Q) = EL=~! P/m^ m (^)- The shape and size of the vesicle 
are determined entirely by the parameters Ro and pi m . The reduced tensors g ab and K a b do 
not depend on R Q . Therefore, Rq can be expressed as a function of A and pi m via the relation 
A — J d£ly/g = RqJ dQy/y. In the disordered phase p — 0, and R = R - We will use the 
Hamiltonian H = TL^ + 7i CU rv expressed in terms of reduced parameters and the constant 
area A in our calculations of shape changes below the second-order disordered-to-n-atic 
transition. 




(1) 




(2) 
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The similarities between H and the Landau- Ginzburg Hamiltonian for a superconductor 
in an external magnetic field, 

H LG = J d 3 x[r\^\ 2 + C|(V - te*A)^\ 2 + ^u\^ + i-(V x A - H) 2 ], (3) 

are striking (Here e* = 2e/hc). Both have a complex order parameter ip with covariant 
derivatives providing a coupling between ip and a "vector potential" A or A a . In a magnetic 
field, the superconductor can undergo a 2nd order mean-field transition from a normal metal 
to the Abrikosov vortex lattice phase with a finite density of vortices determined energetically 
by temperature and the magnetic field H. The magnetic field is conjugate to the vortex 
number N v since / <i 3 x(V x A) = LN v <f> , where L is the length of the sample along H and 
0o = hc/2e is the flux quantum. On a closed surface with n-atic order, there is a second- 
order mean-field transition to a state with vortex number determined by topology rather 
than conjugate external field. Thus, the n-atic transition on a closed surface is analogous 
to transition to an Abrikosov phase with a fixed number of vortices rather than fixed field 
conjugate to vortex number. 

Before proceeding with our analysis of the n-atic transition on a deformable sphere, it is 
useful to recall Abrikosov's calculation of the transition to the vortex state. The first step 
is to calculate the eigenfunctions of the Harmonic part of Hql when V x A = H. These 
can be divided into highly degenerate sets separated by an energy gap huj c = 2Ce*H. In 
the Landau gauge with A = (0, Hx, 0), the eigenfunctions in the lowest energy manifold are 
tpk = e ik y e - e * H ( x ~ x k) 2 w here Xk = k/e*H. The order parameter tp{x) of the ordered state 
is expressed as a linear combination ^(x) = J^Ck^k where the complex parameters Ck are 
determined by minimization of 7Ygl 

We will now proceed to analyze our problem in the same spirit. We first diagonalize Ti,2 
in the high-temperature spherical state when p = and = = — cos 9, A e = Ag = 0, 
i.e., we determine the functions i\) that satisfy CD a D a i\) = eip, for p = A a = A° a . In the 
lowest energy manifold, ip will have exactly 2n zeros specifying the positions of vortices of 
strength 1/n. In the stereographic projection gauge where z = tan(6 ) /2)e i ^ > , we find 
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2M \ n 1 2n 

^=^ttw) v^ z ~ Zi)i (4) 



and e = Cn/R. This function has zeros at the 2n points Zi and only at these points. Since 
these functions are polynomials in z multiplied by a common prefactor, any linear combina- 
tion of them will yield another function of exactly the same form but with different positions 
of zeros. Thus, Eq. is the most general function in the lowest energy manifold. We note 
the similarity between ip 1 / 71 and fractional quantum Hall wavef unctions [11]. Functions in 
higher energy manifolds will have more zeros corresponding to the creation of ± vortex 
pairs. Reexpressing z and 2$ in polar coordinates and choosing a appropriately, we obtain 

ij = tP J\ (sin ~ cos ^ e *(*-*0/ 2 - cos ~ sin h e -*{*-<MI^\ = ^p(^) ( 5 ) 
i=i V 2 2 2 2 J 

as a function of the polar coordinates Qi = 0j) of the positions of its zeros on a sphere. 

Since ip of Eq. (|5|) is the most general function in the lowest energy manifold, the order 

parameter and vesicle shape (determined by p) just below the transition are obtained by 

minimizing 7i over ip and the positions of zeros. To order ifj A , we can write 7i as 



H — I (19. 

4n 



|^[(r-r c )H 2 + iuH 4 ] + Jdn P (n)&(n) + ^Y, l (i 2 -m + 2)\pim\\ (6) 

where r c = —AnnC/A and 



= (A/2n)(r\^ + \u\^) - —M^j") = (7) 

L Sill (7 

with J b = g- 1/2 5H/5A b evaluated at pim = and Rq = R. Minimization over p and -00 
leads to ip 2 = — r[{Qi}]/u[{Qi}}, 

«- = ~ w =m 2) f daY '- im s - ^-f) ( / + 2) *- ( ">- (8) 

and the effective free energy density, 

f\m\]- l r2[m] (9) 

depending only on the positions {fij} of the zeros of P(9) where 
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r[{f!,}] = (r-r c )|df!|F(n)| 2 , (10) 

The next step is to minimize over f2j Though in principle straightforward, this 

task becomes quite complicated as n grows large. In order to evaluate it is necessary 

to integrate |P(fi)| 2 and |P(f2)| 4 over f2 and to calculate p\, m for arbitrary The function 

|P(f2)| 2 can be expressed as Ili n (l — cos7j) where cos7j = cos 9 cos 6>j + sin 6 sin 6>j cos(0 — 
Thus |P(f2)| 2 is a polynomial of order 3 2n in cos 6 and sin#. |$(fi)| is more complicated; 
it is a sum of (n/2)(n + 1) + 1 polynomials of order 3 2n . We have been able to evaluate 
f({Vti}) analytically for n — 1 and n = 2 (where |P(f2)| 2 has 81 terms). For these two 
cases, we find, as expected that the minimum energy configurations are, respectively, those 
with zeros of |P(fi)| at opposite poles and at the vertices of a tetrahedron. For n > 2, a 
complete evaluation of /({fij}) for arbitrary becomes a daunting task. We, therefore, 
appealed to symmetry to treat n = 3, 4 and 6. For n = 3 and n = 6, we expect the zeros of 
|P(Q)| to lie, respectively, at the vertices of an octahedron and an icosahedron. For n = 4, 
following the calculations of Palthy-Muhoray[8], we expect the zeros to lie at the vertices 
of a distorted cube obtained by rotating its top face about its four-fold axis by 7r/4 and 
compressing opposite faces. For this case, we minimized the energy over a single parameter 
determining the separation between opposite rotated faces. 

The shape function p(il) = ^oP*"^ (^) associated with n-atic order is proportional to 
tpQ ~ r — r c to the order of our calculations. In general, the Legendre decomposition of 
will contain Legendre polynomials of order 2n. For n = 1,2, and 3, we find 

?<2,(Si) = llv? |y? + vf ( - i? + 173)1 + - {t ( - y * + Y *~ 3)] 

Fig. 1 shows the shapes described by these functions and the shapes n = 4 and n = 6 
(whose shape functions are too long to display in this letter). The transformations from the 
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initial spherical shape to the distorted shapes occur continuously. Our calculations for the 
shape are valid to order r — r c . As temperature is lowered, higher order terms in r — r c and 
higher order spherical harmonics are needed to describe the equilibrium shape. In Ref. [5], 
a variational function for ip and spherical harmonics up to order 8 were used to calculate 
the shape for n = 1 for temperature well below the transition. The shapes described by the 
functions in Eq. (|Tl|) exhibit unacceptable negative curvature regions at low temperature, 
and higher order spherical harmonics (with the appropriate symmetry) are needed to provide 
a correct description. 

We have presented here an analysis of the mean-field transition to n-atic order on a fixed- 
area surface of genus zero. Our analysis is very similar to that of Abrikosov for the transition 
from a normal metal to a vortex lattice in a type II superconductor at H c2 and completely 
ignores the effects of fluctuations, which may lead to qualitative changes in our results. We 
are really dealing with two kinds of order: n-atic order and the positional order of vortices. 
In mean-field theory, these two kinds of order develop simultaneously. In superconductors, 
fluctuations drive the normal-to-superconducting transition in a field first order [14]. In the 
Abrikosov phase, fluctuations of the vortex lattice destroy[15] superconductivity but not 
long-range periodic order. In two-dimensions, screening of vortices drive the Kosterlitz- 
Thouless transition in an infinite superconductor in zero field to zero[16]. Both the above 
effects may be important for n-atic order on a sphere, and we are currently investigating 
them. 
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FIGURES 



Mean-field shapes of deformable surfaces of genus zero with n-atic order. Above the mean- 
field transition temperature, the equilibrium shape is spherical for all n. Below, the tran- 
sition, the equilibrium shape depends on n and in general has a polyhedral form with 2n 
vertices that coincide with the positions of strength 1/n vortices. 
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